UWO-TH-08/17 



The Fate of the Sound and Diffusion in 
Holographic Magnetic Field 



1,2 



Evgeny I. Buchbinder 1 , Alex Buchel 

1 Perimeter Institute for Theoretical Physics 
Waterloo, Ontario N2J 2W9, Canada 
2 Department of Applied Mathematics 

University of Western Ontario 
London, Ontario N6A 5B7, Canada 

Abstract 

It was shown in [1] that in the presence of the magnetic field the sound waves in 
(2+1) dimensional plasma disappear and are replaces by a diffusive mode. Similarly, 
the shear and charge diffusion fluctuations form a subdiffusive mode. However, since 
the limit of small magnetic field does not commute with the hydrodynamic limit it 
is not obvious whether or not these modes are stable under higher order corrections. 
Using AdS/CFT correspondence we show that in the case of the M2-brane plasma 
these modes do exist as we find the corresponding supergravity solutions. This allowed 
us to compute the conductivity and the shear viscosity to all orders in magnetic field. 
We find that the viscosity to entropy ratio saturates the Kovtun-Son-Starinets bound. 
This extends the universality property of the shear viscosity to the case of the strongly 
coupled plasma in external magnetic field. 
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1 Introduction 



In pQ we studied first-order viscous magneto- hydrodynamics of strongly coupled (2+1)- 
dimensional conformal systems in the framework of gauge theory/string theory corre- 
spondence of Maldacena [2|. The existence of essentially soluble holographic model, 
i.e., the M2-brane plasma [2j [3j 0], allows one to probe intricate aspects of strongly 
coupled relativistic conformal viscous fluids in the presence of external magnetic field. 
In particular, it is known that the hydrodynamic limit in (2+1) dimensions does not 
commute with the limit of small magnetic field [5j. As a result, one expects a drastic 
modification of the transport properties of magnetized fluids. Indeed, in [1] is was 
found that the sound wave in magnetic plasma can propagate only in the limit of 
vanishing magnetic field. Depending on the scaling of the magnetic field in the hydro- 
dynamic limit either only the attenuation or both the attenuation and the speed of the 
sound waves are effected by the background magnetic field. In the latter case, it was 
found that the magnetic field reduces the effective speed of propagating sound modes, 
while enhancing their attenuation. The field theoretical arguments in the setting of 
Hartnoll-Kovtun-Muller-Sachdev (HKMS) magneto-hydrodynamics [5] further suggest 
that a hydrodynamic mode with linear dispersion relation disappears from the spec- 
trum for finite magnetic field. Similarly, in this regime, the standard diffusive modes 
in viscous fluids become subdiffusive, with u oc — iq 4 dispersion. We would like to test 
these predictions in the holographic model of magneto-hydrodynamics of M2-brane 
plasma. 

Our second motivation is to understand how background fields affect the viscosity 
of strongly coupled fluids. Previously, it was discovered that relativistic holographic 
plasma fluids (with various gauge groups, matter content, with or without chemical 
potentials for conserved U(l) charges, with non-commutative spatial directions) have a 
universal value of the shear viscosity at infinite 't Hooft coupling [H [7J 13 EE [101 [11] . The 
universality of the ratio of the shear viscosity to the entropy density extends also to non- 
relativistic holographic CFTs [T2"l IT5] . Additionally, four-dimensional conformal CFTs 
with a dual holographic description and equal central charges a = c have a universal 
leading finite 't Hooft coupling correction [I3l|T5]. O n the contrary, the leading non- 
planar correction to the ratio of shear viscosity to the entropy density is not universal 
[TB] . Since in magnetized (2+1) fluids the shear mode becomes subdiffusive, and the 
dual holographic setting falls outside the most general universality class considered in 
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[9], one naturally questions whether the ratio of shear viscosity to the entropy density 
continues to be universal. 

The paper is organized as follows. In the next section, we review magneto-hydro- 
dynamics of HKMS and its soluble holographic realization as hydrodynamics of dyonic 
black holes in M-theory. In section 3, we extend the supergravity analysis of [1] and 
discuss propagation of 'sound waves' in M2-brane plasma in the hydrodynamic limit 
with finite external magnetic field. In section 4, we study 'shear modes' and compute 
the ratio of the shear viscosity to the entropy density of the M2-brane plasma in the 
external magnetic field. In section 5, we comment on a computation of the shear 
viscosity using the Kubo formula. Some technical details are presented in Appendix 
A. 

After this work was completed, the paper [T7J appeared which approaches magneto- 
hydrodynamics from the gravity side along the lines of [18] . In the future, it would be 
interesting to compare the results of [17] with the HKMS approach. 



2 Magneto- Hydrodynamics and Dyonic Black Hole Geometry 

2.1 Hydrodynamic Modes in the Presence of Magnetic Field 

In this section, we will review magneto- hydrodynamics in (2+1) dimensions follow- 
ing [T]. We are interested in hydrodynamic properties of the (2+1) dimensional theory 
on the large number of M2-branes in the presence of the external magnetic field. This 
theory can be understood as the maximally supersymmetric gauge theory in three 
dimensions at the infrared fixed point. The equations of motion and the conformal 
properties can be rigorously derived using the fact that the theory admits a holo- 
graphic dual description as M-theory on AdS± x S 7 [21 [3l H]. The appropriate field 
theory equations then follow from the symmetries of the AdS background. This was 
studied in detail in [I] and here we will quote the results. The relevant field theory 
equations of motion are just the conservation laws of the form 



f) u T — J? T u 

<V M = o, 



(2.1) 



where T^ v is the stress-energy tensor, J M is the current and F^ v is the external electro- 
magnetic field. In the present paper, it is taken to be magnetic, that is 

F 0i = 0, i = 1,2, I), <,,!'>. (2.2) 
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Another important equation is 

K = > ( 2 - 3 ) 
which means that the field theory under study is conformal. Eq. ( 12.31) is the consequence 
of the fact that the leading near-the-boundary asymptotics of the gauge field in AdS 
is constant, which implies that the magnetic field represents a marginal deformation. 
See [1] for details. 

The expressions for T^ v and J v to first order in derivatives were derived in [5] 
from postulating positivity of the entropy production along the lines of Landau and 
Lifshitz |19j . The stress-energy tensor is given by the standard expression 

T» v = eu»u v + PA"" - riiA^A^idaUp + d p u a ) - A^d^u r ) - (A>* u d a u a . (2.4) 

Here 

A"" = f + «V, (2.5) 

is the fluid 3-velocity, e and P are the energy density and the pressure respectively, 
and rj and ( are the shear and bulk viscosity. Since our theory is conformal it follows 
that £ = 0. It is important to note that P is different from the thermodynamic pressure 

P = p- MB , (2.6) 
where M is the magnetization. Also note that conformal invariance implies that 

Similarly, the current J M is given by 

J» = pu » + a Q A^(-d v fi + F va u a + £d v T) , (2.8) 

where p is the charge density, \i is the chemical potential, T is the temperature and oq 
is the conductivity coefficient. To study fluctuations around the equilibrium state 

vf = (1, 0, 0) , T = const. , fi = const. , (2.9) 

we choose (8u\ = 5u x , 6u 2 = Su y , ST, Spi) as the independent quantities. As usual, 
all the fluctuations are of the plane- wave form exp{—iujt + iqy). In this paper, we are 
interested in hydrodynamics with no net charge density and, correspondingly, with no 
chemical potential 

p = 0, /i = 0. (2.10) 
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In this case, as was shown in [T], the equations for the linear fluctuations get separated 
into the two decoupled pairs. The first pair reads 



w(^) 5T - q(e + P)5u y = 



( dP\ 

u(e + P)5u y - q J ST + iq 2 r]5u y + ia Q B 2 Su y = 



(2.1i; 



If we set B = these equations describe the sound waves with dispersion relation 

uj = ±^--iq 2 —^—, (2.12) 
2 q e + P ' V ; 

where eq. ( 12 .7J) has been used. We will refer to these equations as to the "sound 
channel" . The second decoupled pair of equations is 



u(e + P)5u x — qBoqbp + i<jQB 2 bu x + iq 2 r]5u x = , 
' dp^ 
dp, 



'dp\ (2-13) 
u [ — ) 5p + qaqB5u x + iq 2 Oqbp = . 

T 



Note that here we are assuming that the susceptibility does not vanish at 

p — p — 0. If we set B = these two equations further decouple. One equation 
describes the shear mode 5u x with dispersion relation 

u = -zq 2 -^ J5 . (2.14) 

The other one describes the charge diffusion mode 5p with dispersion relation 

^ = -^ 2 y^-. (2.15) 

We will refer to eqs. ( 12.131) as to the "shear channel" . If we turn on the magnetic field 
B, the hydrodynamic modes undergo a drastic change. The reason, as one can see 
from eqs. (12.111) and (12.131) . is that the limit of small B does not commute with the 
hydrodynamic limit of small u and q. So the magnetic field cannot be thought of as a 
small perturbation. In the case of non-zero B (kept fixed in the hydrodynamic limit) 
we obtain the following solutions. In the sound channel we do not get the sound waves 
anymore. Instead, we obtain a constant solution 

B 2 , v 

<*> = - %G Q-^-p ' (2.16) 
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and a diffusive mode 

iq 2 t + P 

u = ~^m- (2 - 17) 

One can interpret (I2.17P as that the effective speed of sound vanishes once the magnetic 
field is turned on. In the shear channel, the usual shear and diffusive modes disappear. 
Instead, we also obtain a constant solution ( 12. 16ft and a subdiffusive mode 

u = -iq' ^r-. (2.18) 

B 2 ( f£) 

The modes ( I2.17P and (12.181) will be the main focus of our paper. 

Since the limit of small B does not commute with the hydrodynamic limit one can 
worry that the solutions (12.171) and (12.181) cannot be trusted. Indeed, these solutions 
imply a hierarchy of amplitudes. From eqs. (12. lip and (12.171) it follows that 

ST 1 , 
, 2.19 

ou y q 



and from eqs. (12.131) and (12.181) it follows that 

6fi 1 



[2.20) 



Su x q ' 

Hence, given the amplitudes of the linearized fluctuations 8u x and 5u y , we find that 
the amplitudes of the linearized fluctuations 5T and <5/i are strongly enhanced in the 
hydrodynamic limit. Then one can expect that the solutions (12.171) and (12.181) are 
unstable under higher order (higher derivative) corrections to T^ v and J M . That is, 
terms which are naively of higher order because they are suppressed by higher powers 
of oj and q can, in fact, modify hydrodynamic equations at lower order because they 
come with a large amplitude. 

Unfortunately, at the level of the effective field theory it is very difficult to answer 
whether or not these solutions exist. However, for the case of the M2-brane plasma 
we can use the description in terms UV complete M-theory on AdS^ x S 7 background. 
Our gravitational analysis in the later sections will show that in the case of the M2- 
brane plasma with large number of M2-branes the modes (I2.17P and (12.181) do exist. 
Moreover, finding these solutions on the gravity side will allow us to calculate the 
conductivity coefficient <jq and the shear viscosity rj to all orders in magnetic field. 
These results indicate that AdS/CFT correspondence is a helpful method to study 
hydrodynamic modes whose very existence is subtle from the field theory prospective. 
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2.2 Supergravity Magneto-Hydrodynamics 

According to AdS/CFT correspondence, in the limit when the number of M2-branes 
becomes very large, their dynamics can be described by the eleven- dimensional su- 
pergravity on AdSi x S 7 . For our purposes, this theory theory can be consistently 
truncated to Einstein-Maxwell theory on AdS& [20] . The supergravity action is then 
given byE 



S 



— / d 4 £A 



-g 



R _L TP TP MN 
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(2.21) 



where the bulk coupling constant g is given by 



J_ V2N 3 / 2 

~J2 ~ 



gr 6tt 

where N is the number of M2-branes. The corresponding equations of motion are 



(2.22) 



R 



MN 



IFmlF^ — -QmnFlpF lp 



3g 



MN 



V M F 



MN 



(2.23) 







The equilibrium state of magneto-hydrodynamics is described by (asymptotically AdS^) 
dyonic black hole geometry with planar horizon whose Hawking temperature is identi- 
fied with the plasma temperature. The solution looks as follows 



(2.24) 



ds 2 = -d{r) 2 dt 2 + c 2 (r) 2 (dx 2 + dy 2 ) + c 3 (r) 2 dr 2 
F = ha 2 dx A dy + qadr A dt , 



where 



and 



Cl (rf 



a 



f(r), c 2 (r) 2 . 

— i'— T I i y 1 1™ 



c 3 (r) ; 



a 



f(r) = l + (h 2 + q 2 y -(l + h 2 + q 2 y . 



(2.25) 



(2.26) 



In these coordinates, r = 1 corresponds to the horizon and r = is the boundary. The 
black hole parameters (h, q, a) are related to the field theory magnetic field, chemical 
potential and temperature as [21] 



B = ha 2 



-qa 



T 



a 



(2.27) 



1 For simplicity, we set the radius of AdS± to unity. 
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Note that the ^-component of F goes to a constant ha 2 on the boundary and is 
identified with the boundary theory magnetic field B. On the other hand, the t- 
component of the vector potential behaves near the boundary as A t = —qar. It is 
interpreted as the chemical potential in the boundary theory. 

Now we list some thermodynamic properties of the dyonic black hole. See [21] for 
more details. The appropriate thermodynamic potential is obtained by evaluating the 
(renormalized) action ( 12.211) and is given by 



g z 4 V a 1 a 4 



where V is the area of the (x, y)-plane and p is the thermodynamic pressure. The other 
quantities of importance are the density of energy, entropy and electric charge. They 
are given by 

1 a 3 / /i 2 B~ 
g z 2 Y a a 

71 

S = — 
9 

and 



, o + ^ + \ • (2-29) 



2 a 2 , (2.30) 



1 

p = —afx. (2.31) 
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In addition, we introduce magnetization per unite area 

V\dB) Tlt g'a { 1 

Just like on the field theory side, we introduce 

P = p- MB . (2.33) 



One can show [21] that it is P rather than p that coincides with the spatial compo- 
nents (T xx ) and (T yy ) of the stress-energy tensor, just like we have in eq. (j2.4p . It is 
straightforward to check that 

P = \> (2-34) 

which is consistent with conformal invariance. In this paper, we consider magneto- 
hydrodynamics in the absence of the net charge density p. Thus, we set q = 0. Note 
that even though p and \x vanish the derivative 
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is non-zero. 

To study a holographic dual of the hydrodynamic modes, we need to find linear 
fluctuations of the supergravity equations of motion (12. 23f) around the black hole back- 
ground (ESQ), 



QMN QMN + hMN 

Am — * Am + o,m ■ 



(2.36) 



It is convenient to impose the gauge 

In parallel with field theory, the fluctuations gMN and om will be of the form exp(—iut+ 
iqy) and the r-dependence is to be obtained from solving the linearized Einstein and 
Maxwell equations (12.231) . As was explained in [I], for both q and h non-zero, all the 
metric and gauge field fluctuations couple to each other and no decoupling of various 
modes exist. The reason is that the background (12.241) does not have any symmetry 
which usually allows one to decouple scalar-, vector- and tensor-type fluctuations. 
However, in the case of interest q = there exist two sets of decoupled fluctuations. The 
first set corresponds to the field theory sound channel. The corresponding fluctuations 
are 

\htti hty, h xx , hyy, a x } . (2.38) 

The second set corresponds to the field theory shear channel and includes the following 
fluctuations 

{h tX} h xy} a t} a y } . (2.39) 

In the next section, we will show that the fluctuations (I2.38j) indeed correctly describe 
the diffusive mode (I2.17p . In section 4, we will show that the fluctuations (12.391) indeed 
describe the subdiffusive mode (I2.18p . 
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3 The Fate of the Sound Waves 



In this section, we will consider the equations of motion for the fluctuations (12. 38ft . Let 
us introduce 



h a =ci(r) 2 h tt = e~^y cArf H tt 



h ty =c 2 (r) 2 h ty = e 
h xx =c 2 (r) 2 h xx = e 
e 



-iujt+iqy 



c 2 (r) 2 H, 



ty i 



-iujt+iqy 



c 2 (r) 2 H x 



(3.1) 



h yy =c 2 (r) 2 h yy = e-^y c 2 (r) 2 H yy , 



a T =ie 



"yy 

-iujt+iqy 



(I, 



where H tt , H ty , H xx , H yy and, a x are functions of the radial coordinate only and Ci(r) 
and c 2 (r) are defined in eqs. (I2.25P and (I2.26[) . Expanding eqs. (12.231) to linear order 
we obtain the following system of equations 



o =K + H' tt 



111 



qcg 
c 3 



\H XX ~\- Hyyj 



hi 



£2 
Cl 



2c\ 



\ 2 2 

+ u 2 {H xx + H yy ) + 2uq H ty )-3^ h 2 a 4 (H xx + H yy ) + 6 °-\ ha 2 q 

J C 2 C 2 



(3.2) 



=H" y + H' ty 



In 



c 2 



ClC 3 



f 2 2 
Co _ _ C' 



+ -§ uq H xx - 4^ ha 2 [ha 2 H ty + u a x ) 

c 2 c 2 



=H'xx + 2 H 'xx 



In 



qc 2 



1 ff' 



In 



£2 

Cl 



(3.3) 



+ ^2 ( ^ (-^^ ~~ ~~ ? ~2~(Htt + -^kt) 



\ 2 2 

2o;g# t2/ ) - ^ h 2 a 4 (H xx + H yy ) + 2 ^ /ia 2 g 

/ C 2 C 2 



(3.4) 



=iy w + 2 H 'yy 



In 



qc 2 



#1 



In 



C2 
Cl 



+ 



2c 2 



-f^cz) + q 2 -2~{Htt — H x 

C 2 



\ 2 2 

+ 2uqH ty ) - ^| /i V {H xx + H yy ) + 2 °-\ ha 2 q a x 

/ C 2 C 2 



=a x + a4 



In 



£l 

c 3 



a-, a; 



- 2 \ c 2 / 

+ 7T2 h ® 2 l( H tt + H xx + Hyy) 

Z 2 J ^ C 2 \ 



r 2 



+ 2u -y H ty 

c 2 



(3.5) 
(3.6) 
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In addition, we obtain three first class constraints from varying the action with respect 
to the gauge fixed metric components h tr , h yr and h„. 



=uj \H. 



Hyy] + 



ln^ 

Cl 



(H xx + H yy ))+q[ HL + 2 



ln^ 

Cl 



H 



ty > 



+ q — (H tt — H x 



(3.7) 
(3.8) 

(3.9) 

If we set in these equations h = 0, we see that we can also consistently set a x = 0. 
Then the remaining equations for H tt , H ty , H xx and H yy can be shown to coincide with 
those in [22] and describe the sound waves with dispersion relation f)2.12p . See [22] for 
details. To continue, we note that the gauge f!2.37j) does not fully fix the diffeomorphism 
invariance. Since we have five (second order in derivatives) equations and three (first 
order in derivatives) constraints, there are precisely two combinations invariant under 
the residual gauge transformations. They were found in [1] to be 



=q [[H tt - H xx }' - 

c 2 ( 

=[lncic 2 ]' [H xx + Hyy]' - [In 4]' H' tt + -| loo 2 {H xx + H yy ) + 2uq H } 

c i V 

\ 2 2 

) _ h 2 a 4 (H xx + H yy ) + ha 2 q a x . 

J C 2 C 2 



Z H 

Za =a x 



4 i Hfy + 2 // 



yy 



2H V 



q 2 c\c x 



OJ* c 2 c 2 



: 2 C 2 J UJ 2 ci 



(3.10) 



2q 



ha 2 (H xx - H x 



tni i 



Then from eqs. f l3.2p - fl3.6l) and fl3.7p -f l3T9l) we obtain two decoupled gauge invariant 
equations for Zh and Za 



=A H Z" H + B H Z' H + C H Z H + D H Z' A + E H Z A 



(3.11) 



=A A Z" A + B a Z'a + C A Z A + D A Z' H + E A Z H . (3.12) 

The connection coefficients {Ah, • • ■ , Ea} can we computed from eqs. f l3.2p - fl3.9l) and 
fl3.10p using explicit expressions for the q's in fl2.25p . Since these coefficients are very 
long and cumbersome we will not present them in the paper. Below we will present 
these equations in the limit of small uo and q. 

As the next step, we will discuss the boundary conditions. According to the general 
prescription [23l 121] , in order to obtain the dispersion relation (poles in the retarded 
Green's function) we have to impose the following boundary conditions 
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• Zh and Za are incoming waves at the horizon r = 1. 

• Zh and Za satisfy the the Dirichlet boundary conditions on the boundary r = 0. 
That is, both Zu and Za have to vanish at r = 0. 

In [TJ it was shown that Z# and Za have the following behavior at the horizon 

Z H {r) = /(r)-^(r) , 
Z A {r) = f(r)-^ 2 z A (r) , 

where we introduce 

m = £r< ' = 2ST (3 ' 14) 

The functions z# and ^ are now regular and non- vanishing at the horizon. In addition, 
they have to satisfy the Dirichlet boundary conditions at r = 0. 

A crucial ingredient in search for the solution is a proper understanding of the 
correct relative normalization. Since our equations (13. lip and (13.121) are homogeneous, 
we can normalize one of the functions, say zh, to be unity at the horizon r = 1. 
However, after that we cannot normalize za- Moreover, the ratio |— ■ can depend on q. 
Since we are going to solve the equations of motion perturbatively in q it is important 
to establish how — scales with q. We recall that we are looking for a diffusive solution 
with uj ~ q 2 . Then from eq. (13. lUj) it follows that for small q 

z H ~ -zH u , z A ~ a x . (3.15) 

Note that due to an additional symmetry between x and y at q = it follows that 
H xx = H yy at q = and the leading term for small q in za is a x . Furthermore, we know 
that Hu is dual to the tt-component of the boundary stress-energy tensor whereas a x 
is dual to the x-component of the boundary current J x . Going back to the field theory 
side, using eqs. (12. 4p and (12.81) one can show that 

5J X 5u„ 



ST u ST ( 3 - 16 ) 

where eq. (12. 19ft has been used. Then we find that j 4 - ~ g 3 . Let us now parameterize 
the ansatz for our solution. We parameterize the dispersion relation as follows 

xo = -i£c. (3.17) 

The additional h 2 in the denominator is dictated by the field theory result (12.171) . The 
coefficient C is now assumed to have a perturbative expansion in h. Similarly, it is 
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convenient to pull out the appropriate powers in h in the q-expansion of Z H and Za- 
We end up with the following ansatz 



Z H = f{rY™ 12 [f x {t) + f£ F 3 (r) + 
Z A = f(rr^(F 2 (r) + 0(q 2 )) 



(3.18) 



One can show that with these powers of h in the denominators, the functions i^(r) 
have now a perturbative expansion in h 2 . Note that since zh is normalized to be unity 
at the horizon, we can choose Fi(r) = 1 and Fs(r) — at r = 1. 

First, we will solve eqs. ( 13. lip and (13. 12)) to leading order in h, that is ignoring 
the /i-dependence in C and Ft(r). Doing this we will determine that C can be fixed 
entirely by imposing the proper boundary condition on the fluctuations Zjj and Za 
near the horizon. Thus, we will be able to generalize the procedure and find an analytic 
expression for C to all orders in h without actually finding the full analytic solution for 
the fluctuations. We start with the equations to the leading order in q. From eq. (13.111) 
we obtain 

F"--F[ = 0. (3.19) 

r 

The solution with the prescribed above boundary conditions is 

Fi(r)=r 3 . (3.20) 

From eq. (13.211) we get 

2 + r 3 , 3aC 2 

Fo ~ , ox ^2 + = , 3.21 

2 r (i_ r 3) 2 ^8r(l-r 3 ) ' V ; 

where the solution for Fi (13.201) have been used. The general solution to this equation 
is given by 



F 2 {r) = C 1 + C 2 ln(l-r 3 ) 



aC 2 

+ 



32 



(1 + 2r\ 
— j=- J - 2 ln(l - r) + ln(l + r + r 2 ) 



(3.22) 



The integration constant C 2 has to be fixed by requiring that F 2 (r) is regular at the 
horizon. It gives 

C 2 = ^. (3.23) 
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The other integration constant C\ is fixed by requiring that F 2 (r) vanishes at r = to 
be 

„ airC 2 

Note that we are not able to fix the diffusion constant C working at leading order in 
q. We have to go to next-to-leading order in q and consider the equation for 

p „ _ l F , r 4 (81(-4 + r 3 ) - 432C(-2 + r 3 ) + 512C 2 r(-l + r 3 )) = 

3 r 3 9 6 (! _ r 3)2 ' y ■ > 

where the above solutions for F\{r) and F^ir) have been used. It is possible to find the 
general solution for F$(r). We will not write it here because it is rather lengthy. The 
solution has a logarithmic singularity at the horizon of the form 

F 3 (r)~^(-9 + 16C)ln(l-r). (3.26) 

Requiring that the solution is smooth fixes C to be 

C = 1 + 0(h) . (3.27) 

The two integration constants in the solution for F 3 (r) are fixed by requiring that F 3 (r) 
vanishes at r = and r = 1. 

From this procedure it becomes clear that to obtain C we need to understand the 
near horizon structure of the solution and require that it is smooth. It is possible to 
perform this analysis for arbitrary h. The details are presented in Appendix A. This 
allows us to find the exact value of C 

C = ^-(3-h 2 )(l + h 2 ). (3.28) 
lb 

Thus, we have managed to reproduce the diffusive mode in the sound channel on the 
gravity side! We proved that this mode does exist and the exact (to all orders in h) 
value of the diffusion constant is given by eq. (13.281) . 

To summarize the results, we have obtained a supergravity solution with the fol- 
lowing dispersion relation 

«> = -^(3-/> 2 )(l + /> 2 ) + 0(q 4 ). (3.29) 



2 One can worry that going to next-to-leading order in q cannot fix C because the equation for F3 
will also depend on higher order coefficients in the dispersion relation. However, it is straightforward 
to check that it is not the case. 
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Let us compare it with the field theory counterpart (I2.17p . For this we will rewrite (I2.17P 
in terms of (», q,a,h). From eqs. (jSHD , (12T27]) and f[2729|) it follows that fl2TTj) can 
be written as 

iq 2 3 



/i 2 16 



(3-h 2 )(l + h 2 ) 



9 2 °q 



(3.30) 



Comparing (13.301) and (I3.29P we obtain the following answer for the conductivity coef- 
ficient 

1 V2N*/ 2 



(3.31) 



g 2 6tt 

This coincides with the result for oq obtained earlier in [5J [25] to leading order in 
B. Our result (13.311) . however, is valid to all orders in B. Our calculation provides a 
rigorous proof that the conductivity coefficient does not depend on the magnetic field. 



4 The Fate of the Shear Modes and Charge Diffusion 



In this section, we will consider the fluctuations (12.391) describing the shear channel. 
Let us introduce 



h 



tx 



ll X y 6 



-iut+iqy 



c 2 (r) H tx , 

xy ) 



at =i e - iut+iqy a t , 



(4.1) 



a y =ie 



-iut+iqy 



where H tx , H xy , a t and a y are functions of the radial coordinate r. Expanding eqs. (12.231) 
to linear order we obtain the following system of equations 



o =K + 



In 



ClC 3 



H[ x -%{Ah 2 a* + q 2 c 2 2 )H tx - q ^H xy 



r 4 



r 2 



Aha 2 qc\ A Aha 2 ujc\ A 
2 — <k j a 



y ■ 



(4.2) 



In 



ClC 2 

c 3 



i , 2 r 2 



C 



a" + 



cic 3 



g 2 c| 



to; 1 



In 



£1 

C3 



2^2 



-a v + 



at + 



ha 2 ujc\ 



-H, 



tx 



(4.3) 
(4.4) 
(4.5) 
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In addition, we have two first class constraints obtained from varying the action with 
respect to the gauge fixed components h xr and a r 



and 



= A + -£ H l + > ( 4 - 6 ) 



o = 3«; + 3 s * • ( 4 - 7 ) 

If we set /i = we see that the equations describing (H tx , H xy ) and (a t , a y ) decouple 
from each other. In this case, the equations for (H tx , H xy ) describe the shear modes with 
dispersion relation (12.141) [26J. Similarly, the equations for (a t ,a y ) describe diffusion 
with dispersion relation (12.151) [26J. 

To continue, we introduce quasinormal modes invariant under the residual diffeo- 
morphisms 



Z H = qH tx + ujH xy , 

2 (4.8) 
Za = qcit + oja y ha H xy . 

Then from eqs. (j4.2j) - ([4.8p we obtain two decoupled gauge invariant second order equa- 
tions for Zh and Z& of the form (13. lip . (I3.12p . These equations are rather lengthy and 
we will not write them in the paper. Below, we will present them to the lowest orders 
in uo and q. The quasinormal modes Z^ and Z& have the same boundary conditions 
as discussed in the previous section. Namely, Zh and Za are incoming waves at the 
horizon r = 1 and vanish on the boundary r = 0. Repeating the same analysis as in 
the previous section, we arrive at the following ansatz for our solution 



Z H (r) = /(r)--/ 2 ^ (V x (r) + q 2 F 3 (r) + ^F 5 (r) + G(q 6 ) 
Z A (r) = f(r)-^ 2 (F 2 (r) + q 2 F 4 (r) + ^F 6 (r) + 0(q 6 )) 



(4.9) 



where 



and the functions -F,(r) are non-singular and non- vanishing at the horizon and have 
a perturbative expansion in h. In addition, they satisfy the Dirichlet conditions on 
the boundary. Since the equations are homogeneous, we can choose the following 
normalization at the horizon r = 1 

F 2 (r)| r=1 = 1 , F 4 (r)| r=1 = F 6 (r)| r=1 = . . . = . (4.11) 
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Once (14.111) is chosen, no further normalization condition for Fi, F 3 , F 5 , . . . at r = 1 
can be imposed. We parameterize the dispersion relation as follows 

to = -i^C, (4.12) 

where the diffusion constant C is assumed to have a perturbative expansion in h. As 
in the previous section, we will first solve the equations to leading order in h, that is 
ignoring the h dependence in C and i^(r). Then we will generalize our method for 
arbitrary h and find C to all orders in h. To leading order in q we obtain the following 
simple equations 

F!'--F[ = 0, (4.13) 
r 

and 

F 2 " = 0. (4.14) 
The solution with the prescribed above boundary conditions is 

F 1 (r) = bir 3 , F 2 (r) = r . (4.15) 

The integration constant b\ so far is not fixed. Our analysis shows that it is fixed by 
requiring that the function F 5 (r) is non-singular at the horizon. We will not present 
the details of this analysis since they are not important for our purposes. At the next 
order in q we obtain equations for F3 and F4. They can be solved analytically but 
these solutions are not of interest for us. It turns out that the diffusion constant C is 
determined from the equation for F 6 which reads 

pi/ , F „ Cr(27r(-1 + r 3 ) + 16C(2 + r 3 )) 3CV(2 + r 2 ) 

Fe + 9(1^) 2 9(-l+r3)2 F 2 + 2( _ 1+r 3 )2 F2 = • ( 4 - 16 ) 

Note that to determine F 6 we need to know only F 2 which was found in (14.151) . The 
other functions Fi, F 3 , F 4 do not enter this equation. It possible to find the general 
solution for Fq analytically. We will not write it because it is rather lengthy. It turns 
out that to determine C it is enough to study the behavior of Fq near the horizon. We 
find that near the horizon F 6 has a logarithmic singularity of the form 

F 6 (r) = ^(27 - 32C) ln(l - r) + . . . , (4.17) 
54 

where the ellipsis stands for the non-singular terms. Since Fq has to be smooth at the 
horizon it follows that we can have a non-trivial solution 

C = ^ + 0(h). (4.18) 



17 



The two integration constants of _F 6 are fixed by requiring that _F 6 vanishes at r = 
and r = 1. 

Now it is clear how to generalize the above procedure for arbitrary h. The diffusion 
constant C is obtained from requiring that the solution does not have a logarithmic 
singularity at the horizon. In fact, it possible to carry out the near-horizon analysis 
for arbitrary h. Repeating similar steps as in Appendix A, we obtain the final result 

C=<mi. (4.19) 
32 v ; 

Thus, we have managed to reproduce the subdiflusive mode (12.181) in the shear channel 
on the gravity side! We proved that it does exist and found the value of the diffusive 
constant C to all orders in magnetic field in the large N limit. 

To summarize, we have obtained the supergravity solution with the following dis- 
persion relation 

h 2 32 v ; 

Let us compare it with the field theory counterpart (I2.18p . Using eqs. (I4.10p . (12.271) 
and (I2.35P we can rewrite (12.181) in the form 

^SLz^V ( , 21) 

The factor (3 — h 2 ) 3 comes from rewriting the temperature in terms of a and h. Com- 
paring (14.201) and (I4.2ip we conclude that 

" = < 4 - 22 > 



where, to recall, a is related to the temperature and the magnetic field as follows 

m a ( B 2 

T = IT 3 ~ - 
Air \ a 



3 - — ) . (4.23) 



Note that despite the simple form (14.221) . r], after being rewritten in terms of T and B, 
has a non-trivial dependence on the magnetic field. Recalling the expression for the 
entropy density (12.301) we, finally, obtain 

- = 7- ■ ( 4 - 24 ) 
s An 

We see that the ratio rj/s saturates the KSS bound [7]. This extends the universality 
theorem of [9j for the case of strongly coupled plasma in external magnetic field. 
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5 The Kubo Formula 



In the previous section, we extracted the ratio of the shear viscosity to the entropy 
ratio from the dispersion relation of the shear quasinormal modes. Alternatively, one 
can use a Kubo formula, as in j9], to directly evaluate the shear viscosity: 

where G xyxy (G^ xy ) is the advanced (retarded) two-point correlation function of the 
stress-energy tensor (with indicated spatial indices) evaluated at zero momentum. 
Naively, the universality arguments presented in [9] do not hold here. Indeed, the 
dual gravitational mode, h xy does not generically decouple as one can see from eqs. 
f l4.2l) -( ]4~5]) . Furthermore, the background of the bulk vector field (holographically dual 
to the background boundary magnetic field) is not solely polarized along the time di- 
rection, as assumed in [9]. In this section, we revisit the argument of [9] and extend 
the universality theorem to the case of background magnetic field in (2+l)-dimensional 
strongly coupled plasma. 

First, notice that even though the gravitational mode h xy does not decouple for non- 
zero momentum q, see eqs. (I4.2l) - ()4.5p . the decoupling occurs for q = 0, which is all 
what is needed for the computation of the correlation functions in (15.11) . Physically, the 
reason why such a decoupling occurs is because for vanishing q there is an additional 
symmetry in plasma associated with a reflection along the y-axis. As a result, the 
graviton polarization h xy is the only fluctuating mode which is doubly parity odd 
under reflections along the x- and ?/-axis. Hence, it must decouple. Second, from (14.31) 
we see that the rescaled (see eq. (14.11) ) graviton wavefunction H xy at q = satisfied the 
equation of motion for the minimally coupled massless scalar in a regular Schwarzschild 
horizon geometry (I2.24j) . But now, we are precisely in the setup of the universality 
arguments of [9] ! One can literally repeat the analysis presented there to establish that 
the shear viscosity rj, as determined by ( 15.11) . is proportional to the entropy density s 
with the proportionality coefficient as in (I4.24H . 
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A Exact value of C 

For a background magnetic field h held fixed in the hydrodynamic limit, the sound 
channel dispersion relation is given by ( 13. 17ft . with the constant C = C(h). Here, we 
outline steps necessary to obtain exact analytical expression for C(h). 

Within the ansatz (13181) . F x is given by (13T20D . From (IBTTTj) and (13TT2D we further 
find the following equations for F 2 and F 3 

= F"+ 2 ^h 2 -r\h 2 + l)-2 pi _ ( h2 - ?>)*C 2 ; (A1) 

2 r(r 3 h 2 — r 2 — r — l)(r — 1) 2 8r(r 3 h 2 — r 2 — r — l)(r — 1) ' 

r 4 /i 2 -l 12r 2 /i 2 „ 

=F" - 6—-— — -F' + 



r(r 4 /i 2 -3) d r 4 /i 2 -3 a 
r 2 (/i 2 - 3)(r 3 (/i 2 + 1) - 4)(r 4 /i 4 + r 4 /i 2 - 16h 2 r + 9h 2 + 9) 

2C 2 a(r 3 h 2 - r 2 - r - l)(r - l)(r 4 /i 2 - 3) 2 ( 

24r 2 (/i 2 - 3)/i 2 



C 2 a(r 4 fr 2 - 3) 



F 2 + Jc 



3 • 



where 



* -32 (r -l) W -X-r-l)W-3) ( 512r3<r " 1)(rV " ^ " '' " 1)C2 
+ 16r 2 (-54(/i 2 + 1) + 84/i 2 r + 27(/i 2 + l) 2 r 3 - 60h 2 (h 2 + l)r 4 + 24/i 4 r 5 

+ 3/i 2 (/i 2 + 1) V - Qh 4 (h 2 + l)r 8 + Ar 9 h 6 )C - (h 2 - 3) 2 (-36(/i 2 + 1) + A0h 2 r 

+ 9(h 2 + 1)V - 32h 2 (h 2 + l)r 4 + A8h A r 5 + h 2 (h 2 + 1)V - \2h\h 2 + l)r 8 

+ 8r 9 /i 6 )^ . 

(A.3) 



It is straightforward to construct power series solutions first for F 2 ( IA.1I) and then for 



^3 ( j A. 21) near the horizon. Regularity of {i<2, F3} for small x = 1 — r then uniquely 



fixes C as in (I3T28D . 
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